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Abstract We consider a measurable stationary Gaussian stochastic process. A criterion for 
testing hypotheses about the covariance function of such a process using estimates for its norm 
in the space L p ( T), p > 1, is constructed. 
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1 Introduction 

We construct a criterion for testing the hypothesis that the covariance function of 
measurable real-valued stationary Gaussian stochastic process X(i) equals p{r). We 
shall use the correlogram 


T 



0 < t < T, 


0 


as an estimator of the function p(r). 
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A lot of papers so far have been dedicated to estimation of covariance function 
with given accuracy in the uniform metric, in particular, the papers [2, 4, 6, 11, 
12] and the book [13]. We also note that similar estimates of Gaussian stochastic 
processes were obtained in books [7] and [1], The main properties of the correl- 
ograms of stationary Gaussian stochastic processes were studied by Buldygin and 
Kozachenko [3], 

The definition of a square Gaussian random vector was introduced by Kozachenko 
and Moklyachuk [10]. Applications of the theory of square Gaussian random vari¬ 
ables and stochastic processes in mathematical statistics were considered in the paper 
[9] and in the book [3]. In the papers [5] and [8], Kozachenko and Fedoryanich con¬ 
structed a criterion for testing hypotheses about the covariance function of a Gaussian 
stationary process with given accuracy and reliability in ^(T). 

Our goal is to estimate the covariance function p(r) of a Gaussian stochastic 
process with given accuracy and reliability in L p ( T), p > 1. Also, we obtain the 
estimate for the norm of square Gaussian stochastic processes in the space L p ( T). 
We use this estimate for constructing a criterion for testing hypotheses about the 
covariance function of a Gaussian stochastic process. 

The article is organized as follows. In Section 2, we give necessary information 
about the square Gaussian random variables. In Section 3, we obtain an estimate 
for the norm of square Gaussian stochastic processes in the space L p ( T). In Sec¬ 
tion 4, we propose a criterion for testing a hypothesis about the covariance function 
of a stationary Gaussian stochastic process based on the estimate obtained in Sec¬ 
tion 3. 


2 Some information about the square Gaussian random variables 
and processes 

Definition 1 ([3]). Let T be a parametric set, and let 5 = {£ t : t £ T} be a family of 
Gaussian random variables such that EG = 0. The space SGs(O) is called a space 
of square Gaussian random variables if any ( £ SGs(O) can be represented as 

C = G r Ae-E^ T AG 

where £ = (£ 1 , • • •, Gv) T with £& £ ■H', k = 1,..., n, and A is an arbitrary matrix 
with real-valued entries, or if £ £ SG~(f2) has the representation 

C= lim (^Af n ~E^AQ. 

n—f 00 v 7 

Theorem 1 ([3]). Assume that £ £ SGs(f2) andYaiC, > 0. Then the following 
inequality holds for \s\ < 1; 

Eexp {7lte)} s iA^l exp {^} = tl,w (1) 

Definition 2 ([3]). A stochastic process Y is called a square Gaussian stochastic 
process if for each t £ T, the random variable Y (f) belongs to the space SG~(f2). 
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3 An estimate for the L P (T) norm of a square Gaussian stochastic process 


In the following theorem, we obtain an estimate for the norm of square Gaussian 
stochastic processes in the space L p (T). We shall use this result for constructing a 
criterion for testing hypotheses about the covariance function of a Gaussian stochastic 
process. 

Theorem 2. Let {T, 21, / i}be a measurable space, where T is a parametric set, and 
let Y = {Y(t),t £ T } be a square Gaussian stochastic process. Suppose that Y is 
a measurable process. Further, let the Lebesgue integral J T (Ey 2 (f)) ? dp(t) be well 
defined for p > 1. Then the integral J T (Y(t)) p dp(t ) exists with probability 1, and 


P\ / I Y(t)\ P dp(t) > £ > < 2 


'\ 


t , £l/p V / 2 

1 H-t— exp<j 

pp 
Gp 


£P 


V2 Cl 


( 2 ) 


for all e > + y'd + 1 )p) p C p , where C p = f T (EY 2 (t))? dp(t). 

Proof. Since max I> o x a e~ x — a a e~ a , we have x a e~ x < a a e~ a . 

If C is a random variable from the space SGs(O) and x = , where 

s > 0, then 


E 


f-LJiLV 


< a°e~ a ■ Eexp 


\V2fE( 2j 


and 


Eicr < 


(f2EC 


a a e~ a Eexp 


V s 

From inequality (1) for 0 < s < 1 we get that 

s C 


/-?_£L\ 

\V2fEC) 


Eicr < 


(y/ZE? 


< 


V s 

2 ( a/2EC 2 


a a e~ a Eexp 


VI — s \ s 
= 2L 0 (a)| 


a a e “exp< — 


V2 yfEC- 

s 

7 ! 


+ Eexp< — 


5 C 
V2 v/E? 


YV2E7 


V 5 


(3) 


Let Y(t), t £ T, be a measurable square Gaussian stochastic process. Using the 
Chebyshev inequality, we derive that, for all l > 1, 


P 


J\ Ydp(t) > £| 

T 


E(J T \Y{t)\ p dp(t)Y 

£ l 


Then from the generalized Minkowski inequality together with inequality (3) for 
l > 1 we obtain that 
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< J (2L 0 (s)(2T]Y 2 (t)) 2 (pl) pl s pl exp{—pl}) 1 dp(t) 

T 

= (2L 0{s) yJ (2EF 2 (f)) = s~ p {pl) p exp{-p}dp(t) 

T 

= (2L 0 (s)) J 2% s~ p (pl) p exp{-p} J (E Y 2 (t))*d/j,(t). 

T 

Assuming that C p = / T (EF 2 (i))^ dp(t), we deduce that 

E|F(i)| p d/z(f)^ < 2L 0 (s)2%(lp) pl exp{- P l}C l p s~ pl . 

T 


Hence, 


p[J | Y(t)\ p dp(t) > ej < 2- (2iyL 0 (s)(p p ) l (exp{~p}) l C l p (s- p ) 


u^i ( l p Y 


= 2 L 0 {s)a l (l p )\ 

V 1 — 

2 2" C • — 2 ^ 25(7 p 

where o = „ p , that is, a? = — i -T £ -. Let us find the minimum of the function 

ePsPe ’ ’ -±- 


t/j(7) = a l (l p ) 1 . We can easily check that it reaches its minimum at the point l* = 

l 

eoP 

Then 


-V / 1 \f—r -V 

2Lo(s)il>(l*) = 2Lo(s)ae°- p ■ | —— ) eaP = 2Lo( s ) aeaP 

\ea p J 

= 2L 0 (s)expj— PCSl ' L 1 =2L 0 (s)exp|- 

^ 2JpeC p J L 

'1 e 1 ! p 


S£ p 

2 5 cl 




ex Pi -s 7t 


2^Cp 


In turn, minimizing the function 9(s) = 2 exp{—s(P H—r)} i n we 

V S 2 ic p 


1 + 


^ttf- Thus, 

^1/P 


6>(s*) = 2 


\ 


1 + 


cj 


exp< - 


£P 


v^Cj 


deduce s* = 1 — 
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Since l* > 1, it follows that inequality (2) holds if —V > 1. Substituting 

ea P V2pC p 

the value of s* into this expression, we obtain the inequality £ 2 ' p > pC 1 p /p (C 1 p /p + 
\/2e 1//p ). Solving this inequality with respect to e > 0, we deduce that inequality (2) 
holds for e > (-^= + + 1 )p) p C p . The theorem is proved. □ 


4 The construction of a criterion for testing hypotheses about the covariance 
function of a stationary Gaussian stochastic process 


Consider a measurable stationary Gaussian stochastic process X defined for all t £ R. 
Without any loss of generality, we can assume that X = {X(t), t £ T = [0, T + A], 
0 < T < oo, 0 < A < oo} and E X(t) = 0. The covariance function p(r) = 
E X(t + r)X(t) of this process is defined for any ret and is an even function. Let 
p(r) be continuous on T. 

Theorem 3. Let the correlogram 

T 

p( t ) = J X(t + T)X(t)dt, 0 < T < A, (4) 

o 

he an estimator of the covariance function p{r). Then the following inequality holds 
for all e > (-4L + y/(^ + l)p) p C p : 


j (/SO) - p{r)) P dT > ej < 2 
o \ 


e l /Py/2 

1 4- 1 — exp 


a 


p 



where C p = J^(T—u)(p 2 (u) + p(u + T)p(u — T))du)^dTandO<A< oo. 

Remark 1. Since the sample paths of the process X(t) are continuous with proba¬ 
bility one on the set T, p(r) is a Riemann integral. 


Proof. Consider 

E(p(r) - p(r)) 2 = E(/5 (t)) 2 - p 2 {r). 

From the Isserlis equality for jointly Gaussian random variables it follows that 


T T 

E(/5(r)) 2 -p 2 (r)=E^ J j X(t + T)X(t)X(s + T)X(s)dtd S ^-p 2 (T) 

0 0 


1 


T T 


E X(t + r)X(t)EX(s + t)X(s) + E X{t + r)X(s + r) 


o o 


x EX{t)X(s) + E X(t + r)X(s)EX(s + r)X(t))dtds - p 2 (r) 
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1 

J^2 


T T 


0 0 
T T 


p 2 {r) + p 2 (t — s) + p{t — s + r)p(t — s — r))dtds — p 2 (r) 


= T 1 I I <yf>2<yt ~ S )+P( t - S + T )P T )) dtds 
2 

— J^2 


0 0 
T 


f(T- u)(p 2 {u) + p{u + r)p(u — r))du. 


We have obtained that 


1 

E (/S0) - p{r )) 2 = ^2 J (T ~ u)(p 2 (u) + p(u + r)p(u - t) ) du. (5) 
o 

Since p{r) — p(r) is a square Gaussian stochastic process (see Lemma 3.1, Chapter 6 
in [3]), it follows from Theorem 2 that 


p{t) - p{T)) P dT > £ \ < 2 


■\ 




■ exp 


cip 


£P 




Applying Eq. (5), we get 


C p = 


A T 

2 


T 2 / 
o o 


{T — u)(p 2 (u) + p(u + r)p{u — r))du\ dr. 


The theorem is proved. 
Denote 


□ 


ff(e) = 2 


\ 


1 + 


z 1 /pV2 

cl 


■ exp 


£P 


V2 Cl 


From Theorem 3 it follows that if £ > z p = C p (-^= + y/{% + l)p) p , then 

A 

p { / ( /5( - r ' ) ~ p ^y dr > e } - 9 ^>- 
0 

Let £s be a solution of the equation g(j£) = 5, 0 < 6 < 1. Put Ss = maxje^, z p }. It 
is obvious that g(Ss) < S and 


A 

p { / _ p ^y dT > - 5 - 


( 6 ) 
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Let H be the hypothesis that the covariance function of a measurable real-valued 
stationary Gaussian stochastic process X(t) equals p(r) for 0 < r < A. From The¬ 
orem 3 and (6) it follows that to test the hypothesis H, we can use the following 
criterion. 

Criterion 1. For a given level of confidence 5 the hypothesis 11 is accepted if 

A 

J ( P ( T ) ~ p( T )) P dp(r) < Sg; 
o 

otherwise, the hypothesis is rejected. 

Remark 2. The equation g{e) = S has a solution for any S > 0 since g(e) is a de¬ 
creasing function. We can find the solution of the equation using numerical methods. 

Remark 3. We can easily see that Criterion 1 can be used if C p —> 0 as T —> oo. 
The next theorem contain assumptions under which C p —► 0 as T -A oo. 

Theorem 4. Let p{r) be the covariance function of a centered stationary random 
process. Let p(r) be a continuous function. If p(T) —>■ 0 as T — > oo, then C p —> 0 as 
T —» oo, where C p = J^(fi(T,T)) p ^ 2 dt and 

T 

T ) = f ( T ~ u ) + p ( u + T ) p ( u ~ r ))du, A > 0, T > 0. 

o 

Proof. We have ip(T, r) < tp f^(p 2 (u) + p(u + r)p(u — r))du < 4p 2 (0). Now it 
is suffices to prove that i/j(T, r) —> 0 as T —> oo. From the L’Hopital’s rule it follows 
that 


T 

lim ib(T,T) = lim — / (p 2 (u) + p(u + r)p(u — r))du 
T-too T^ooT J 

0 

= r lim (^(T) + p(T + r)p{T - r)) = 0. 

Application of Lebesgue’s dominated convergence theorem completes the proof. 

□ 

Here are examples in which we find the estimates for C p . 

Example 1. Let H be the hypothesis that the covariance function of a centered mea¬ 
surable stationary Gaussian stochastic process equals p(r) = B exp{— o|t|}, where 
B > 0 and a > 0. 

To test the hypothesis H, we can use Criterion 1 by selecting [>t(t) that is defined 
in (4) as an estimator of the function p{r). Let 0 < A < oo. We shall find the value 
of the expression 

T 

0 
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i i i 

= jre-^du + T f e-^e-^'du - J ue'^du 
0 0 0 


- J ue~ a \ u+T \e~ a \ u ~ T ' [ du 

o 


= I\ + I2 + ^3 + I4 ■ 
Now lets us calculate the summands: 


h = 


1 

T [ e- 2au du=^(l-e~ 2aT ), 
J ^ 


Io=T e" a| “ +r| e" a| “" r| ^ 




= T[ I e~ 2aT du + / e~ 2au du 


= t( re~ 2aT - — e~ 2aT + — e~ 2aT 
\ 2 a 2 a 

T T 

h = [ ue~ 2au du = —-^—e~ 2aT + 2. [ e~ 2au du 
J 2a 2a J 

0 0 

= _ZLp- 2aT _L p -2aT , J_ 

2a 4a 2 4a 2 ’ 


L= ue- a| “ +r| e- a| “- T| ^ 


r _/ 

J ue- a(u+T) e a(u - T '>du + J ue- aiu+T) e~ a(u - T) di 


— 2ar , 


^ 2a “du 


r 2 J' 

= —e~ 2aT - —e~ 2aT 
2 2a 


_Lp-2ar _ 1 —2aT , 1 „-2aT 

2a 4a 2 + 4a 2 
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Therefore, 


I =[Tt + 


T 


<[Tt + 


2 a 
T 
2 a 


t 2 t 1 

T ~~ 2^ ~~ la 2 


-2ar | _J_ -2aT . ^ ^ 

2 a 2 2 a 4 a 2 


-2ar i ZL + J_ e - 2aT 

+ 2a + 2a 2 


Thus, we obtain 




2a 


T 


1 


Tt+ — e~ + — + — 


2 a 2 a 2 


p/2 


dr 


P T P//2 P 1 

= m^^rh = ( 2Bp^I 5 , 


where J 5 =/ 0 A ((t+ i)e 2aT + £ + 2aT )P/ 2 dr. 

Example 2. Let H be the hypothesis that the covariance function of a centered mea¬ 
surable stationary Gaussian stochastic process equals p(r) = B exp{—a|T| 2 }, where 
B > 0 and a > 0. 

Similarly as in the previous example, to test the hypothesis H, we can use Cri¬ 
terion 1 by selecting pt{t) debited in (4) as the estimator of the function p(r). Let 
0 < A < oo. Let us hnd the value of the expression 


1 

1 = J(T- u)(e~ 2au2 + e- a|u+r| V a|u - T|2 )du 

0 

T T T 

= j Te~ 2au2 du + T J e - a|u+r| V a|u - T|2 du-y ue~ 2au2 du 


- J zte -a| “ +T|2 e -a| ‘ u-r|2 dzi 

o 

= h +12 + h + Ia- 

Now let us calculate the summands: 


h=T [ e~ 2au2 du < T [ e~ 2au2 du = ^ T 


2 V2a’ 


h=T e" o|,i+r| V a|u " T|2 du = Te" 2ar2 / e" 2a “ 2 du < 


V^T _ 

2\J 2 a 


2 CLT Z 


1 1 

h = J = -j-J e-^4-^) = -j- a (e- 2 ‘ T ‘ ~ 1 ), 
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T 


T 



,u = e 


— 2 ar 2 


/ 


ue~ 2au2 du = ^l e - 2a " 2 (e- 2aT2 - 1). 

An v ' 


0 


0 


Hence, 



Thus, we obtain 


V 


A 


)) V '\r= m 


Cp< 



p 1 



0 
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